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Abstract 

We study a one-component quaternionic wave equation which is 
relativistically covariant. Bi- linear forms include a conserved 4-current 
and an antisymmetric second rank tensor. Waves propagate within the 
light-cone and there is a conserved quantity which looks like helicity. 
The principle of superposition is retained in a slightly altered man- 
ner. External potentials can be introduced in a way that allows for 
gauge invariance. There are some results for scattering theory and for 
two-particle wavefunctions as well as the beginnings of second quan- 
tization. However, we are unable to find a suitable Lagrangian or an 
energy-momentum tensor. 
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1 Introduction 



Many attempts have been made to consider the extension of the usual quan- 
tum theory, based upon the field of complex numbers, to quaternions. The 
1936 paper by Birkhoff and von Neumann pQ opened the door to this pos- 
sibility, and the 1995 book by Adler [2] covers many aspects that have been 
studied. 

Here is a wave equation that appears to have escaped previous recognition. 

dip 

— i = u • vip + rmjjj. (1.1) 

The single wavefunction ip is a function of the spacetime coordinates x, t. 
The usual elementary quaternions i, j, k, are defined by 

i 2 = j 2 = k 2 = ijk = -1 (1.2) 

and 

„ d d , d , 
u-V = i—+j— + k—. 1.3 

ox ay oz 

Boldface type is used to designate a 3-vector. 

This combination ()1.3|) of elementary quaternions and space derivatives 
was orginated by Hamilton [3] in 1846; it's square is the negative of the 
Laplacian operator. 

What one should note about the equation (jl.lj) is that it employs quater- 
nions which multiply the wavefunction on both the right side and the left side. 
This distinction arises from the non-commutativity of quaternion algebra and 
is central to the present study. 



2 Other Equations 

There are other quaternionic wave equations one can consider, based upon 
the apparent structural similarities between quaternions and relativity. The 
simplest is 

which, when squared, appears as a 4-dimensional Laplace equation, not a 
wave equation. 
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Going to two-dimensions we construct 




= m*. (2.2) 

When this equation is squared, we do get a wave equation, but it is for a 
tachyon. If one sets m = in this equation, it can be revised to appear as 
either two copies of the Weyl equation or the Maxwell equations (keeping 
only the imaginary components). 

Various authors have shown that the familiar Dirac equation can be put 
into quaternionic form. This may be done by putting an i to the right of 
\1> on one side of equation (|2.2|) 4J or by the use of bi-quaternions in a 
one-component equation All of those representations involve eight real 
functions - as does the usual Dirac equation - while the basic equation of the 
current study (jl.lj) involves only four real functions. 

There are two other known relativistic equations with four real compo- 
nents. One of these is the Majorana representation of the Dirac equation, 

iYd^ip = mip (2.3) 

where all four of the gamma matrices can be made purely imaginary, so 
that one can take all four components of the Dirac wavefunction to be real 
functions of space-time. Indeed, if we write our quaternionic wavefunction 

as 

if) = V>o + + 3^2 + kfa, (2.4) 

and arrange these four real functions as a column vector, then our equation 
can be put in exactly this Majorana-Dirac form. One awkward feature 
of that formalism is that the usual Dirac Lagrangian becomes useless for an 
action principle, since every single term is identically zero. 

The other comparison involves the Weyl equation (two complex compo- 
nents), which is usually reserved for massless particles. One can map quater- 
nions onto a two-dimensional space of complex numbers. The correspondence 
can be expressed in terms of the familiar Pauli matrices u — > —ia; and the 
wave equation can be written in a pseudo-Weyl form as 

i— = -ia -W + ma 2 *'. (2.5) 
at 

In this second example, one also has trouble with the usual Lagrangian in 
that the mass term is identically zero. 
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Both of these equations, Majorana-Dirac and modified Weyl, are used in 
the building of supersymmetry theories (see, for example, $]), but only after 
one introduces a second set of wavefunctions - with "dotted" spinor indices. 
Thus, they do return to eight real functions, which are, furthermore, not 
simply real functions but elements of a Grassmann algebra. 

These comparisons leave me without a definitive answer to the question of 
whether the focal equation of this paper (Jl.lj) is truly something new in theo- 
retical physics. The work presented here will be to explore this quaternionic 
wave equation on its own terms and see what interesting things arise. 

3 Some Properties 

In the usual quantum mechanics there is "gauge invariance of the first kind" : 
we can replace the complex wavefunction ip by exp(i6)ip. This freedom is 
also noted by saying that there is a ray, not just one vector, in Hilbert space 
corresponding to each physical state. (The reader will note that this paper 
focuses entirely on the wavefunction approach to quantum theory and not the 
Hilbert space version.) For the quaternionic wavefunction we have a larger 
set of freedoms: ip — > qiipq2, where the two numbers qi,q 2 are quaternions 
of unit magnitude. The one on the left induces a change of basis in the 
elementary quaternions u seen in the equation (jl.ljl . while the one on the 
right changes the particular choice of i and j acting to the right of ip in 
that equation. Thus, instead of the usual U{1) group, we appear to have 



A first calculation is to take another time-derivative of equation and 
arrive at the second-order wave equation, 



which is the ordinary Klein-Gordon equation for a relativistic particle of mass 
m. 

Now we look at some bi-linear forms. The first is p = ip*ip, where the 
complex conjugation operator (*) changes the sign of each imaginary quater- 
nion (and requires the reversal of order in multiplication of any expression 
upon which it operates). The second is the vector U = ip*\xip. While p is 
real, U is purely imaginary and we can write U = iUi +.7U2 + k\J 3 , in terms 
of three real three-vectors. 



SU(2) x SU(2)/Z 2 . 



V 2 ip — m 2 ^ 




dt 2 
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Making use of the wave equation ([1.1)1 . we then calculate 



| = V. Ut , (3.2) 

which is the familiar statement of a conserved current. We shall return to U2 
and U 3 shortly. (If you ask what singled out Ui as the conserved current, it 
is the choice of the imaginary i sitting beside the time-derivative in the wave 
equation (jl.l)) .) 



4 Space-time Symmetries 

Now we look at the behavior of the wave equation (jl.lj) under familiar sym- 
metry transformations. To achieve rotation of the spatial coordinates x, we 
make the transformation 

ip — ► e R il), R = u-6/2 (4.1) 

where 6 is the axis and the angle of rotation. 

For the Lorentz transformation, we start with the infinitesimal form, 

ip — >ip + Bipi, 5 = u-v/2 (4.2) 

where v is the direction and amount of the velocity boost. Note the ap- 
pearance of the imaginary i acting on the right of ip in this transformation. 
I leave it as an exercise for the reader to show that this transformation of 
ip does indeed induce the familiar Lorentz transformation of the spacetime 
coordinates in the wave equation (jl.l)) . 

One can now readily show that the components of the conserved current 
(p and Ui) tansform as a Lorentz 4- vector. With a bit more work, one can 
also see that the other two vectors U2 and U3 transform as the components 
of an antisymmetric second rank tensor in 4-dimensions (also called a six- 
vector). 

A useful notation for operators that may multiply quaternionic functions 
on the right or on the left is the following. 1 

(a\\b)ip = a b, (a\\b) (c\\d) = (ac\\db), (4.3) 

1 A similar notation was introduced by the authors of reference 
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which allows us to write the finite Lorentz transformation operator as e^ B ^ . 
The generators of the Lorentz Group may be constructed as: 

1 d 1 

j = xxV u, K = x— + tV u i. (4.4) 

2 dt 2 11 v 1 

One can extend this to the full Poincare group by adding the displacement 
operators: <9 M = (<9t, V). 

In Appendix A is a more extensive study of various tensors that can be 
built from solutions of the wave equation. 



5 More Bilinear Forms 

Start by defining the derivative operator which acts in both directions, = 
(do, d) = \(d ^ — d M ). This is a covariant four- vector, but let us now see 
how things behave when we combine it with the Lorentz transformation of 
the wave function: 

= i)*d^ -> Dp + - {i, v ■ u^} . (5.1) 

The expression D^ is pure imaginary; and so we can write D^ = Di ^ i + 
^2,^ j +^3,/i k. I hope that the use of the subscripts (1,2,3), denoting which 
imaginary component they come from, does not cause confusion with the 
vector or tensor subscripts /i. The expression inside the anti-commutator 
brackets (next to i) is real. This leads us to conclude that, under the Lorentz 
transformation of ip: 

D 2 ^ and D 3 fl are unchanged (5-2) 
Dx tfl — > Di tfl + more complicated stuf f. (5.3) 

This means that under the full Lorentz transformation of both coordinates 
and wavefunction D 2 and D 3 behave simply as 4- vectors. The quantity Di, 
however, will be shown in Appendix A to be part of a higher rank tensor. 

Before proceeding, we note that D^ =Q can be reexpressed by using the 
wave equation (jl.lj) : 

£> M =o = - ■ U 3 + k[mp + ~V • U 2 ], (5.4) 
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where r = ip*u ■ dip is a real 3-scalar. Under the Lorentz transformation of 
the wavefunction, we calculate r — > r + v • Di. 
We have the identity 

d^Dy = 0; (5.5) 

and we will be interested in the following time derivatives, which are derived 
by using the wave equation (jl.lj) . 

|r— V.D, (5.6) 
= i[2mD 2itx - V • (if)*ud^)} - 2mjD^ + [i, ■ d^] (5.7) 

d 

— U = i[Vp + 2^*u xdip + 2mU 2 ] 

+j[-2mUi + 2D 3 - V x U 3 ] + fc[-2D 2 + V x U 2 ]. (5.8) 
See Appendix A for a more systematic discussion of tensor quantities. 



6 Plane Waves 

One way of representing "plane-wave" solutions of the wave equation (jl.lj) is 

■0(x, t) = exp(r]u ■ p p • x) (f) exp((ir]p + km)t) (6-1) 

where 77 = ±1. The set of possible momentum vectors p = pp should cover 
only one-half of space, to avoid overcounting if solutions. With this, one can 
construct the solution for the general initial value problem: 

rf3x ' Ylj H J^y exp(r]u-pp- (x-x')) 

il}(x',t' = 0) exp((irjp + km)t) (6.2) 

where the subscript H reminds us that the integral covers only half of mo- 
mentum space. 

With the expansions 

exp(rju ■ p p • (x — x')) = cos(p • (x — x')) +r)U ■ p sin(p ■ (x — x')) (6.3) 
exp((ir]p + km)t) = cos(ut) + (ir/p + km) sm(ut)/uj, (6.4) 
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where u = \Jp 2 + m 2 , we sum over rj and reduce (|fi.2|) to the following. 




d 3 p 



2[cos(p • (x - x')) i/j (x', 0) (cos(ut) 



h (2tt) 3 



+kmsm(u>t)/uj) + u • p sin(p ■ (x — x')) ^(x', 0) i sm(ut)/u)]. (6.5) 

Here we can recognize that the results of the integrals over p (which now 
may be extended to cover the full momentum space) give us functions of the 
invariant R 2 = t 2 — (x — x') 2 , which vanish outside the light-cone (R 2 > 0). 
Thus we do have relativistic causality for this quaternionic wave equation; 
something which we could have expected because the solutions satisfy the 
Klein-Gordon equation. 

The Klein-Gordon equation also has the property that positive (negative) 
frequency solutions propagate only to positive (negative) frequency solutions. 
For the quaternionic equation, we have no way to talk about this distinction 
between positive and negative frequencies; however we do find a substitute 
"selection rule" for wave propagation here. 

First, we note the orthogonality relation, 



where I have not required that both sets of momentum variables belong 
to the same half-space. Next, we use this orthogonality in equation ([6.2)1 . 
where we represent ip(x.', 0) as any superposition of plane wave solutions with 
exclusively rj' = +1 (or exclusively -1). The resulting ?/>(x, t) will contain only 
that same value for 77. It is tempting to call this "helicity conservation" in 
the propagation of these quaternionic waves. 

This interpretation is bolstered by the following observations. The op- 
erator u • V, acting on a plane wave solution (|6.1|) . has eigenvalue — r\p. 
Furthermore, one can readily show, from ([l.ip . that 



7 Superposition 

In the usual (complex) quantum theory, if we have two solutions to the 
Schrodinger (time-dependent) equation, ipi and ip2, then any linear combin- 
mation Ciipi + c-^i is also a solution, for arbitrary complex numbers c\ and 





(6.7) 
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c 2 . With our quaternionic wave equation (jl.lj) . the idea of superposition 
requires a slightly different wording. 

Note that the general plane wave solution (j6.1j) has an arbitrary amplitude 
(f) positioned in the midst of certain quaternionic functions of space and time. 
Given any such solution, we find another solution by changing the amplitude: 
<f> — ► q(f>qt, where q and q! are arbitrary quaternionic numbers. Furthermore, 
if we have one solution of (jl.lj) - if)\ with amplitude <j>\ - and another solution 
- ip2 with amplitude 02 - then we also have a solution by simply adding these 
two: if>i + ip2- This version of the principle of superposition is implicit in 
equation (j6.2j) . 



8 Adding Potentials 

The original wave equation (jl.lj) can be extended by the introduction of 
external potentials, as follows. 

dip 

-—% = u ■ X7ip + e</70 — eu • Aipi + mipe ieW j (8-1) 

where ip, A, W are real functions of spacetime. The gauge transformation 
that leaves this equation invariant is: 

<P 
A 
W ■ 

One can show that the previously discussed symmetries still hold, with 
(ip, A) a Lorentz 4-vector and W a scalar. This appearance of the 4-vector 
potentials is (almost) exactly like the usual way of introducing electromag- 
netism into quantum theory; however, the explicit appearance of a gauge 
quantity W is something different. 

The reflection symmetries of equation (j8.1j) are: 

if) — > ip j, t, A, W change sign (T) (8.6) 
if) — ► if> k, x,<f,W change sign (CP) (8.7) 
if) — > if) i, t, x, A, ip change sign (TCP). (8-8) 



if) — > ipe 



■ex 



dx 

A 

W 



2X- 



(8.2) 

(8.3) 

(8.4) 
(8.5) 
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The current conservation equation (J3.2)) is still true for this extended wave 
eqation (|8.1|) ; however, equation (|6.7|) must be modified. For the situation 
where W = 0, we calculate 

-1 = ^-(ip*n ■ dip) = -V • D x + epV • A - 2eA • 0*u x dip) (8.9) 

where we have used the notations from Section 5. From this we see that in 
the case where the only external potential is tp, then the space integral of r, 
which we identified with helicity, is conserved. 



9 More on Plane Waves 

The plane wave solutions to the wave equation (jl.lj) . which we set out in 
Section 6, contain an amplitude which we should study some more. 

■0(x, t) = exp(r) u • p p • x) exp(Qt) (9.1) 

where Q = Quj = (irjp + km) , u = y/p 2 + m 2 . 

We can ask to evaluate the various bilinear forms discussed earlier in the 
case of this plane wave solution. The easiest are 

p = 0*0, t = -VPP, A) = (9.2) 

but to do more we must be able to evaluate 0*u • ]50. 

I now propose to classify the constants in a particular way. The set <p a 
are defined such that they perform a specific rotation, as follows: 

p ■ u a = (f> a Cl (9.3) 

which sends one unit imaginary quaternion into another. With this type, the 
solution can be written as 

ip = 0a exp(Cl(ujt + r/p ■ x)) (9.4) 

which looks like the sort of plane waves we are used to. It should be noted 
that this definition of <p a is not unique but leaves us with a U(l) class of 
equivalent amplitudes, 

0a - ► 0a exp(9Q), (9.5) 
just as in ordinary (complex) quantum theory. 
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With this a type of amplitude, we can now evaluate the plane wave values 
for the following bilinears: 

U • p = iprj — h kp — , D = ?7p(U-p). (9.6) 

Components of the vector U which are orthogonal to p will oscillate rapidly 
in space; thus any space average of them will be vanishing small. 

Two other categories for the amplitudes 0, called (3 and 7, can be defined 

as 

U • P 0/3 = 0/3 j (9.7) 
u ' V 07 — 07 3&- (9-8) 

Note that the three numbers Q, j ,jtt are mutually anticommuting quater- 
nions. If we calculate any of the bilinears involving u • p with either the (5 or 
7 type of amplitude, the result will be rapidly oscillating in time; thus any 
time average will be vanishing small. 

If we stay with the a type amplitudes, we get the following values, in 
the plane wave states, for various 4-vectors that are defined in Section 5 or 
Appendix A. 



j li = p(l,rip/u) (9.9) 

V, = p( V p,u;p) (9.10) 

£> 2 , M = (9.11) 

D^ = pm(l, V p/Lu) (9.12) 



The 4-vectors j and D 3 look like what we would expect for the usual energy- 
momentum. The 4-vector V, however, is spacelike, not timelike; it is similar 
to the spin vector — £^,v,k,\P u S K in the usual theories, where Sq is the 
helicity. 

The plane wave solutions are characterized by a parameter p which we 
sometimes call "momentum". This is merely a linguistic habit carried over 
from conventional quantum theory (following deBroglie's rule that momen- 
tum equals Planck's constant divided by wavelength) and should not be 
confused with the physical quantity called momentum until and unless that 
connection is established. 
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10 Scattering 



Lets start with the wave equation plus a source 

d 

j-ipi = u ■ Vip + mipj + s(x, t)i, (10.1) 
and write the retarded solution as 

^(x,t) = C dt' [ d 3 x' [ -^^y2exp(r]u-pp-(x-x.'))s(^',t')exp(n(t-t')). 

J-oo J JH {27T) d 

(10.2) 

For the general scattering problem, we replace the source s with Vtp and 
add in the initial (free particle) solution ?/>o(x, t). If the interaction V is 
independent of time, then we have an integral equation, 

V>(x, t) = V>o(x, t) + J ^ dt' J d 3 x' ex P(v u • PP • (x - x')) 

V(x.')tlj(x.',t')exp(n(t - If)). (10.3) 

Now we make the "Born approximation" that ip — ipo under the integral 
and let the time t go to +oo. Then, we find that the integral over t' gives 
us a 5{uj — ujq), which is usually read as conservation of energy. This result 
appears to be generally true, not just in the first Born approximation. One 
can now project this solution onto any plane wave solution and achieve the 
quaternionic version of the S-matrix. 

In the special case when the scattering potential V comes from the term tp 
in the extended wave equation (J8.1|) . we also find - as a result of the integral 
over t' - that we have the selection rule n = i] . This is consistent with the 
result noted after equation (|8.9jl . 



11 Some Other Solutions 

We can write solutions for the extended wave equation (J8.1|) in some special 
cases. 

One may ask whether there is a central potential, ip(r), which leads to 
bound states. The easiest way to explore this is through "reverse engineer- 
ing": write down a plausible wavefunction and see what potential fits the 
wave equation. The form 

i;= [/(r)+u-f g(r)](f>(t) (11.1) 
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leads to the requirements 



4/ 2 = 4( r V), ep = -f/9, 0(t)=0 o e fcmi . (11.2) 
dr dr 

If we try the asymptotic (r — > oo) behavior, g — > ar~P, we find a similar be- 
havior for /, provided that < (3 < 2. The wavefunction is then normalizable 



for (3 > 1.5; and the potential is e<p(r) = J (3(2 — (3)/r at large r. Looking 
instead at r — > 0, one can do the same analysis and require (3 < 1.5; this 
suggests we are dealing with something like a shielded Coulomb potential. 

There are familiar procedures for taking the non-relativistic limit of the 
Klein Gordon or Dirac equations. Here is the best I could do with the present 



relativistic equation. First, write i/j = i/j nr exp(kujt) } where u = \Jm 2 + p 2 
m — V 2 /2m. Next, multiply the equation from the right with iexp(—kujt). 
Finally, drop all terms that oscillate rapidly in time, as exp(±2kut). The 
resulting version of the full extended equation (|8.1j) is 

^knH nr ip nr (11.3) 
H nr = -V 2 /2m + m(l - cos(eW)) -en- A\\k (11.4) 

which looks like an ordinary Schrodinger equation except that the single 
imaginary is called k instead of i; and there is also the unfamiliar term with 
A. What looks like an effective potential energy term (coming from the gauge 
quantity W) is positive, thus incapable of producing bound states, although 
it might conceivably yield metastable states through the oscillation of the 
cosine function. 



12 Two Particle Equation 

Previous studies of quaternionic quantum theory have gotten into trouble 
when they try to write wavefunctions for multi-particle systems. In the or- 
dinary (complex) theory, one simply makes a direct product of one-particle 
wavefunctions; and because all the numbers there commute, one can manip- 
ulate such a product to achieve various sensible results. In the quaternion 
case that approach leads to a horrid mess. (See, for example, Adler's book, 
Chapter 9.) 
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The present work suggests a somewhat different approach. Consider this 
construction with plane waves: 

ip(l, 2) = expfau ■ pxp a • xi) ip(2) exp^iU) (12.1) 
ip(2) = exp(r} 2 u ■ P2P2 ■ x 2) <P exp(tt 2 t 2 ) (12.2) 

which might be described as a "nested" product. The symbol here is a 
quaternionic constant, which can depend on all the parameters of this two- 
particle wavefunction. Note that we have written this with independent time 
variables for the two particles. 

These two-particle wavef unctions, with all their momentum-helicity la- 
bels, form a complete orthogonal set of functions in the space of xi and X2. 
Note, however, that this product is ordered in a way that was meaningless 
in ordinary (complex) quantum theory but requires some extra bookkeeping 
in the quaternionic case. 

Let's introduce some more compact notation for such wavef unctions. 

-0 op (l) = exp(r]iu • piPx • xi) 1 1 exp(Qiti) (12-3) 

where the || symbol separates those things that are to act on the left from 
what is to act on the right of whatever follows. Then the two-particle wave- 
function (|12.1|) can be written simply as 

^(1,2) = ^(1)^(2)0; (12.4) 

and we can also write the operator of the wave equation as 

d 

V = — + u- V||i-mp. (12.5) 

Next we have the propagators. 
G(x — x', t — t') = exp(r]u ■ pp ■ (x — x'))| \exp((irjp + km){t — t')) (12.6) 

where J2p,rj = Ih (f^fs an d this leads to 

PG(x - x', t - t') = 0, G(x-x',0) =5 3 (x-x') (12.7) 
G + (x - x') = 6{t - t')G(x - x, t - t') (12.8) 
VG+(x-x') =5\x-x'). (12.9) 
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The coordinate x stands for the full space-time coordinates t, x. Now the 
equation (jl(J.2|) can be briefly written as 



Following that construction, we now write down a general two-particle 
quaternionic wavefunction, as follows 



ip(xx,X2)= / d 4 x[G + (xi — x[) / d 4: x' 2 G + (x2 — x' 2 )s(x' 1 , x' 2 ). (12.11) 



This is a two-particle wave equation of the Bethe-Salpeter type, involving 
separate times as well as separate space coordinates. The term s might be 
left as an external source or might be used to represent some interaction, such 
as V(l, 2)^(1, 2). Note that the order in which the two differential operators 
are applied is significant. 

It seems easy now to extend this to any number of particles. This ap- 
pears to be a significant advance over previous studies of quaternionic wave 
equations, although there are still many issues to be faced. 

13 No Lagrangian 

If I use the interacting wave equation ()8.1j) . and think that ip* is something 
independent of ip, then the following would be suggested as a Lagrangian 
density. 

C = — iip*u ■ Vip — itp*eiptp + itp*eu ■ Aipi — imifj*if)e ieW j. (13.1) 

Varying iip* gives immediately the full wave equation for ip. Before varying 
ipi on the right, we do a few things: partially integrate in space and time; 
move i from left side to right side in the second and third terms and rearrange 
the i and j coefficients in the last term; this is justified because those ip* . . . ip 
expressions are real. Then we get the adjoint wave equation. 

But that prescription is not what the usual action principle allows. The 
familiar game from complex qm does not work here. If one varies each of 




(12.10) 



Acting on this with two of those differential operators gives 



V 2 Vi ip(xi,x 2 ) = s(x 1 ,x 2 ). 



(12.12) 
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the four real functions which make up both quaternionic functions ip and 
ip*, then we actually get 12 equations from the action principle. This is due 
to the fact that this Lagrangian is imaginary, that is, it consists of three 
imaginary parts and each of those parts must vanish after the variation. If 
we write 

£ = i£ 1 +j£ 2 + k£ 3 , (13.2) 

we find that the first term, C%, is Lorentz invariant (see equation (|A.9p ): but 
what should we do with the other two terms? 

Our difficulty with a Lagrangian is different from the difficulty noted 
earlier for the Dirac-Majorana equation or for the pseudo-Weyl equation. 
But we do have a problem here. 

14 Discussion 

Several advances have been made in trying to develop a sensible quantum 
theory based upon quaternions, rather than complex numbers. So far, this 
work has been limited to the wave equation formalism. 

We have noted the lack of a conserved energy-momentum tensor (see 
Appendix A) as well as the lack of a Lagrangian. Nevertheless, we can write 
down the time-development operator as 

U(t)=e m , H = -u-V\\i + m\\k. (14.1) 

This operator H commutes with the angular momentum operator J, but 
whether we want to call it the Hamiltonian is unclear. Perhaps these ques- 
tions wait for a full model of how this quaternionic wave system interacts 
with other physical systems. 

Another approach that may be relevant to that problem, as well as to 
improving our treatment of many-particle systems, is the method of second 
quantization. We are led to write down a quaternionic quantum field operator 
as 

?/>(x, t) =22 ex P{v u • p p • x) a P)T] exp((ir]p + km)t) (14.2) 
p,n 

involving some kind of annihilation/creation operators a p ^ v . With this we 
immediately get 

N = t^t- 3 J d 3 xp = Y: al v a P!V (14.3) 
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k = (2tt)3 j d3xT = -Y, a l,v a P>v VP- ( 14 - 4 ) 

Can one be sure that the matrix product a)a is real? If these are matrices 
in a Fock space of the sort we are familiar with, with nonzero elements only 
on one line parallel to the central diagonal, then this product is real. 

It remains unclear to this author whether the equation studied in this 
paper is merely an alternative mathematical formulation of things already 
well known or whether it may have consequential applications to some as- yet 
unidentified physics. 
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Appendix A: General Tensors 

We can construct Lorentz covariant tensors of any rank, as follows. Start 
with the direct product of the "two-way" derivative operators: 

^ = ^i^ 2 ---^ (A- 1 ) 

where the subscript \i now stands for the set of indices \i\ . . .// n . This ex- 
pression is manifestly a covariant tensor of rank n as far as the coordinate 
transformations are concerned; and our task is to package these between the 
wavef unctions, which transform under infinitesimal Lorentz transformations 
as given in (|4.2|) . 

We note that the packages tp* d^ip are real for n even and imaginary for 
n odd; while these conditions are reversed when we add the quaternions u 
inside the package. 

We find the following constructions for covariant tensors of rank n+1, 

Q { ^ 1] = Q^ = (Q^Q,y- 

Q^ = (VX^{y.V'*^uV}) for n even (A.2) 

= for n odd. (A.3) 

In the case n = this is just the 4-vector current, previously written as 
= (p, Ui). All these tensors are real. 
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In addition, for n odd, we have the tensors of rank n, 

= {y , Rs,» = {y , rd^}, (A.4) 

which generalize the previously noted 4-vectors D 2 ,D S . And for n even, we 
have the tensors of rank n+2, S^x = —S^\ M : 

J 

S^,o,a = S„, a ,o = {—,ip*d^u a il;}, S^p = e aAl {— ,0*^11^} (A.5) 

where a, [3, 7 = 1,2,3. This generalizes the previously noted six- vector 
(U 2 ,U 3 ). 

We can make lower rank tensors by contracting indices: 

g^Q(n+l) = _ (m 2 + I 9 A 9A)Q (n~l) (A _ g) 

where /ii and /12 are in the set /x; and the set fjf has these two indices removed. 
An alternative is to contract one of the /i indices with the v index. We find, 
for solutions of the free wave equation (jl.lj) : 

g^Qin+i) = -mRf-V f or n even (A.7) 

If we are looking to find a Lorentz scalar, a tensor of rank zero, take a 
closer look at the second rank tensor Q^ v . The contraction is 

= {y, rdoif)} + ■ = D lfi + r (A.9) 

and this is exactly zero for the free equation (jl.l|) but not for the extended 
equation ()8.1|) . where it equals —ej^A^. 

Now we look at the contraction of such tensors with the derivative opera- 
tor. In what follows we shall limit ourselves to solutions of the free equation 
(|1.1|) . It is transparent that <9 m Q^ u = for any \i\ in the set of labels //. 
The same holds true for the tensors R and S. Furthermore, by using the wave 
equation, one can show that 

d u Q^u = for n even (A. 10) 

d u Q^ u = 2mRQ for n odd. (A. 11) 
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For tensors of rank one, we have just the previously identified j^, -D 2 ,/i,and 
D 3 ^, all of which are conserved. 

At rank two, the usual desire is for a conserved symmetric tensor, which 
one can call the energy-momentum tensor. The closest we come here is the 
which is not symmetric and is conserved only on the first index. Nev- 
ertheless, this does allow us to write integral quantities which are conserved 
(their time derivatives vanish), as follows. 

V V = (V ,V) = Jd 3 xQ$ (A. 12) 

V = - J d 3 x ij*u ■ dip (A.13) 

V = J d 3 x (Di + mU 3 - x Ui) (A.14) 

This is not what we would identify as the energy-momentum, as noted at the 
end of Section 9. 

For the second rank antisymmetric tensor we have 

S fl ,v = -2mj v + 2D 3>u (A. 15) 

= W,a (A.16) 
d»S^ = -2D 2>v . (A. 17) 
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